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ABSTRACT 

A technique is described for simulating the random modulation im- 
posed by atmospheric scintillation and transmitter pointing jitter on 
a direct detection optical communication system. The system is cap- 
able of providing signal fading statistics which obey log normal, beta, 
Rayleigh, Ricean or chi-squared density functions. Experimental 
tests of the performance of the Channel Simulator are presented. 
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CHANNEL SIMULATION FOR DIRECT DETECTION 
OPTICAL COMMUNICATION SYSTEMS 
M. Tycz and M. W. Fitzmaurice 


INTRODUCTION 

The rapidly increasing need for wideband communication systems for space-to-space and 
space-to-ground links in the late 1970’$ and early 1980’s has stimulated an increased interest 
in optical communication technology. There are. active NASA and DoD supported programs 
which are directed toward developing high data rate, nominally 300 MBps to 1 GBps, optical 
communication systems for space flight use. A typical system might employ a Nd:YAG mode- 
locked laser transmitter using on-off-key (OOK), binary polarization (BPM) or quaternary 
modulation formats and a photoemissive direct detection receiver incorporating either high 
speed, gated photomultiplier tubes or avalanche photodiode/low noise, wideband amplifier 
combinations. Such systems are governed by Poisson statistics [ 1-2] . Their performance for 
OOK or BPM modulation formats and deterministic average receiver signal intensity have 
previously been analyzed and reported in the literature [3-4] . 

Optical communication systems look most attractive in the space-to-space link. In this con- 
figuration, it is easiest to take full advantage of their relatively small size and weight, their 
rugged and long-lifetime components and their low prime power requirements. However, it 
is in this configuration, for example low-earth orbiter to synchronous satellite link, that the 
operational system may require sub-arcsecond pointing of the laser transmitter beam. As a 
result, a dominant noise source may be the interaction of the transmitter far field irradiance 
profile and the instantaneous transmitter pointing direction (see Figure 1 ). Thus, the average 
receiver signal level may be random and obey beta statistics [ 5-7] . For the space-to-ground 
link one can take advantage of the transmitter low power, light weight, etc. system parameters, 
but must contend with the perturbing effects of the atmosphere (see Figure 1). Performance 
analysis reported in the literature for the space-to-ground link incorporate the effects of 
atmospheric scintillation assuming a log-normal channel [8] . Although system performance 
analyses have not been completed for the Ricean and Rayleigh channels, theory does exist 
asserting that certain communication links will be subject to atmospheric scintillation result- 
ing in time varing average signal amplitudes described by Ricean or Rayleigh probability 
distributions [9]. 

Coupling the advancement of the theoretical analyses of direct detection optical communica- 
tion systems for varied channels with the rapidly progressing systems development, it now be- 
comes necessary to design a laboratory device for controlled testing of the communication 
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Figure 1 . Direct detection optical communication channels 


system. It is required that the unit simulate the channel effects of the different link con- 
figurations. This publication addresses the design and performance of one such Channel 
Simulator. 


CHANNEL SIMULATOR SYSTEM CONCEPT 

The purpose of the Channel Simulator is to provide the controlled effects of atmospheric 
scintillation and transmitter pointing inaccuracy for testing and evaluating direct detection 
optical communication systems for the space-to-space and space-to-ground links. As seen 
in Figure 2a, the Channel Simulator transforms a constant signal I into a preselected random 
varying signal I M . The probability distribution /j (i M ) governing the modulated signal level 
of the Channel Simulator output I M (see Figure 2b) will depend on the link to be simulated. 
The Channel Simulator described in this publication is designed to produce random signals 
whose time varying irradiance may obey log normal, chi-squared or beta statistics and whose 
time varying amplitude may obey log normal, Ricean or Rayleigh statistics. It consists of two 
major components: the linearized optical modulator and the processing electronics. Its 
transfer function is given by I M = C V where I M is the irradiance of the modulated light beam, 
C is a constant determined by the modulator crystal parameters and V is the voltage applied 
to the modulator. In general, V may be deterministic or random. If V is random with prob- 
ability density f v (u), then the probability density governing I M will be [10] 


. . fy <W C ) 

j C |. 


( 1 ) 


OPTICAL SUBSYSTEM 

The optical subsystem consists either of a linearized acoustooptic modulator or a linearized 
electrooptic modulator. One such system utilized by the experimenters was a KDP electro- 
optic modulator. Employed in the transverse mode of operation, the modulation field was 
applied normal to the direction of light propagation. Using the axes defined in Figure 3, the 
input polarizer was in the X'Z plane at an angle of 45° from the Z-axis. The output polarizer 
was rotated 90° relative to the input. The transfer function for this modulator system is 
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Figure 3. Etectrooptic modulation - transverse mode. 
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where to is the frequency of incident light, c is the speed of light, n Q and n e are the ordinary 
and extraordinary indices of refraction, E z is the applied field, r 6 3 is the electrooptic coeffic- 
ient of the KDP crystal and (n Q - n e ) L y » oj/2c, the retardation due to natural birefringence, 
is a function of temperature and can be compensated. Therefore, consider 
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If we define E z = V/D and let 
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where is the voltage required from maximum transmission through the crystal. 
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If we perform a series expansion of the right hand side of equation 5, then 
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It is obvious that for V < V 1/2 

^TRANS ^ V2 

— = — T- (7) 
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One can now form a linearized optical modulator system by simply inserting a square root 
module between the applied voltage V and the modulator crystal, assuming operation in the 
range V < V I/2 . . ... ... . * ..... . .. • 

The linearized aeoustooptic modulator is designed using the same prihcipal. Modulation of 
hght due to the interaction of a light Wave with a sound wave induces a diffracted light wave 
at the sum or difference frequency. The transfer function for a beam incident into a medium 
of a distance £ is given by [12]* 



( 8 ) 


where X Q is the wavelength of the incident wave in a vacuum, p is the photoelastic constant 
of the medium, I AC0USXIC the acoustic intensity, u s is the velocity of sound in the acoustic 
medium and p the mass density. Defining I Q as a reference acoustic intensity yielding maxi- 
mum deflection 
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and noting that y/l is proportional to V the applied voltage, then 
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Using the same expansion as before and assuming voltages V < V Q , then the transfer function 
becomes linear where once again a square rooting device is inserted before the modulator 
crystal 
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Results of this technique applied to the Channel Simulator will be presented later in this 
publication. 

PROCESSING ELECTRONICS 

The analog processing electronics generate the random voltages necessary to drive the 
optical modulator for simulating the desired channels. As seen in Figure 4a, it consists of 
Gaussian random signal generators fed into antilog converters, squaring modules and similar 
devices which perform appropriate mathematical operations. Although certain operations are 
shared by more than one channel, the transfer function for each channel will be presented 
separately for clarity. 

Atmospheric Log Normal Channel 

For the atmospheric scintillation channel log normal mode, the random variable that obeys 
log normal statistics is I/I 0 or the instantaneous irradiance divided by the mean irradiance. 

To simulate the log normal behavior of this channel the following random variable I M need be 
generated: 


i M = I o e2L < 12 > 

where 2L is Gaussian with mean K and variance X 2 [13]. The transfer function of the 
channel simulator log normal channel (see Figure 4b) from the Gaussian input, monitored at 
test point TP7, to the modulated irradiance output is given by 

bC -2(X TP7 + A) In 10/ 10 

I M = 10 e 

where A is the log normal offset, b is the gain established by the last amplifier stage, X TP7 
is a zero mean Gaussian random voltage with variance , and C is the linearized modu- 

lator subsystem constant of proportionality. 

Redefining variables in equation 1 3 results in equation 1 2 where I Q = b C/ 10 and L = 

- (X Tp? + A) lnl0/10. Additionally, when simulating the log normal channel, the log 
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Figure 4 , Channel simulator processing electronics (a) composite of all channels (b) log normal atmospheric channel (c) Raleigh or 

Ricean atmospheric channel (d) beta free space channel. 
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Figure 4 (continued) 































amplitude must have a mean equal to negative of its variance s 2 
E {L} = - Var {L> = - s 2 


[ 14] . Therefore, letting 
we establish relationships for the values of A and Uy . with 


A Vs 2 10 /In 10': . 

(14) 

= s 2 [10/ln 10] 2 . 

TP7 

(15) 


The value of s 2 is selected from atmospheric scintillation data and ranges nominally from 
s 2 = 0.01 to s 2 = 0.1 [15]. 


Atmospheric Rayleigh and Ricean Channel 

For the atmospheric scintillation channel Rayleigh or Ricean mode,. the random variable that 
is described by these statistics is the light amplitude. Until now, however, the linear transfer 
function of the Channel Simulator has been referenced to the irradiance level of the output. 
The light amplitude transfer function is thus ; 

;v. 

A=[CV]* = C a V* (16) 


where C, = C 1/2 . 

a . 

A Rayleigh or Ricean random variable P can be generated by the following 

P = [X 2 + X 2 ]^ (17) 

where for the Ricean case X, and X 2 are independent Gaussian random variables with means 
/ij and ju 2 and variances a 2 = a 2 - For the Rayleigh case, ju 1 =p 2 = 0 [16]. 

The total transfer function of the Ricean* channel from the Gaussian inputs to the. modu- 
lated light amplitude output is given by (see Figure 4c) - 

; = C a V A (18) 

where Xj-p, and X rP5 are zero mean independent Gaussian random variables with variance 
<j 2 , and d is the Ricean offset. Note also that for Rayleigh simulation V and therefore the 
modulated irradiance is a chi-squared random variable with two degrees of freedom. 


♦Since it is more general, we present this channel with the understanding that for the Rayleigh case the off- 
set is simply set tozero. . .. • = 


P * 


(Xtp 2 + d) Xppj" 
+- 


20 


20 
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Free Space Channel 

For satellite to satellite optical communications, signal fades caused by pointing jitter may 
limit receiver performance. In this case, the modulation on the average received signal level 
will be described by beta statistics [5-7] . In simulating this fading condition the following 
must be satisfied: 


r M ~ J o e 


-p 2 /2r 


(19) 


where I Q corresponds to the peak of the spatial Gaussian irradiance profile, % is the transmitter 
beam half angle at the 1/e 2 points and p is the distance the beam is from the center of the 
receiving aperture. For typical systems, p 2 is a chi-squared random variable of two degrees 
of freedom corresponding to two independent and orthogonal Gaussian tracking error 
signals. 

The total transfer function of the Channel Simulator in the beta mode of operation can be 
seen in Figure 4d. From the Gaussian random voltage inputs to the modulated irradiance 
output the relationship is 

I M = 10Ce“* ,n 10 / 5 H(Xy P2 + x^,p S )/20] (20) 

where X TP2 and are zero mean, equal variance, Gaussian random variables. Redefin- 

ing I 0 = 10C 

a 2 - Vai^Xppjj. = Var-^Xpp^ 

| 2 = 50/ln 10 (21) 


and 


= 


x| P2 + x^ P5 ) 


then equation 1 9 is satisfied and the density of I M becomes 

t 2 / , J 2 I° 2 


fl M 0m) 


■ 5 © 


o„> { lc 


1 


0<i M <J 0 . 


( 22 ) 


The ratio of £/o is determined by both system design and the link to be considered and will 
vary typically from 3 to 10 [7] . 


COMPONENT TRANSFER FUNCTION TEST RESULTS 

Tests have been performed on the Channel Simulator processing electronics and optical 
modulator subsystem to verify the accuracy of the mathematical operations, the linearity of 
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the optical modulator and the independence of the Gaussian sources. Figure 5 shows the re- 
sults of tests performed on the electronic modules which perform the required mathematical 
operation within the processing electronics. The equations representing the resjiitant curves 
were derived from the data assuming minimum mean square eiror fit to a straight line, least 
square fit to a power curve and least squares fit to an exponential curve where appropriate. 

The transfer function of the linearized acoustooptic modulator system used by the authors 
was tested using a helium-neon laser input beam and an S-20 photomultiplier detector. Fig- 
ure 6 shows typical test results when the modulator bias and the photomultiplier tube voltage 
are set for linear operation. Since the slope of this transfer function is directly affected by 
the photomultiplier tube gain as well as the modulation index, this test was performed be- 
fore each of the channel simulations. 

As previously stated, operation of the Channel Simulator for the beta, Rayleigh and Ricean 
modes assumes that the Gaussian random voltage inputs at test points TP2 and TP5 are zero- 
mean, equal variance and independent. Offset adjustments as well as variable gain stages are 
provided to ensure that the signals are zero mean and equal variance. Tests for independence 
were performed using the Hewlett Packard Correlator, Model 3729A. Recall that two random 
signals X and Y are independent if their crosscorrelation coefficient p x Y is zero, i.e., 


e{xy} -e{x} e{y} 

p = — - — = o 

VVar{x} Var{Yj 


where for this case, the means E{x} and E{y} are zero and the VAR{x} - VAR{y} is 
simply the autocorrelation of X or Y evaluated at a time delay of zero, and E{XY} is the 
crosscorrelation evaluated at a time delay of zero. The auto and crosscorrelation of Xj, p 2 
and X TP5 were taken for Ricean offset d = 0. Using these values at delay zero and substi- 
tuting into equation 23 resulted with the crosscorrelation coefficient p = 0.006. 


STATISTICAL FUNCTION TEST RESULTS 

The usefulness of the Channel Simulator is measured in terms of its ability to accurately 
simulate the desired channel. Figure 7 demonstrates the performance of the log normal 
channel simulation as compared to experimentally measured channel perturbations. The 
plot of log amplitude versus its cumulative probability is on a log normal cumulative proba- 
bility graph. The solid line represents typical results from a GEOS B satellite experi- 
ment [15] for which an argon laser beacon was transmitted through the atmosphere 
to a low earth orbiting satellite. The instantaneous signal was detected on-board the 
spacecraft by a photomultiplier tube and was telemetried to ground for statistical analyses. 
The dashed lines represent data taken from the Channel Simulator system in which the signal 
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Figure 5* Transfer functions of mathematical operations. 




Figure 5 (continued) 
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Figure 6, Transfer function of the linearized acoustooptic modulator. 




Figure 7. Comparison of channel simulator log-normal statistics with GEOS-B data. 

level was log normally modulated. The linearity of the data indicates the modulated signal is 
indeed log normally distributed and the slope specifies the strength of atmospheric scintilla- 
tion. The latter can be made to vary from a log amplitude variance of 0.01 to 0.08. 

Performance of the beta channel requires operation of all the mathematical functions in the 
processing electronics. The output signal of the Rayleigh and Ricean channel is the input 
signal provided to the antilog converter for generating a beta distributed random variable. 

Thus accurate simulation of the beta distribution necessitates good performance of the Ricean 
and Rayleigh atmospheric channel modes. 

Recall that the probability density function governing the beta mode of operation was given 
by equation 22. Therefore, the cumulative density will satisfy 
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Taking logarithm of both sides gives 


2 2 - 

V 1 "' i logl M''T lo8l o <2S) 

Plotting the logarithm of the cumulative beta density versus the logarithm of the relative 
modulated signal level irradiance results in a linear curve, the slope of which is the parameter 
of the beta distribution to be simulated. Figure 8 shows a typical plot of simulated beta 
modulation for which the ratio |/cr = 5.3. The Channel Simulator can be made to simulate 
the space-to-space link with beta perturbations ranging nominally from 1 < £/cr < 1 0. 

Although not directly applicable to simulating optical communication channel effects, an 
additional mode of operation has been coupled to the system. By adding an optional square 
rooting stage to the output of the Rayleigh and Ricean channel, the voltage output and, 
therefore, the modulated irradiance (as opposed to the modulated E-field) becomes Rayleigh 
or Ricean distributed. A test was performed on the operation of this Rayleigh random voltage 
output. Recall that a Rayleigh distributed random variable P has cumulative density 

F p(P>=^-/* 27rre" r2 / 2<j2 dr) > 0 (26) 

2310 Vo 

where a 2 is the parameter of the Rayleigh distribution. Performing the integration in equa- 
tion 26 results with 


V- 1 n [1 - F (p)] 



(27) 


Figure 9 is a plot of equation 27 using the square root output voltage for data. The linearity 
of the data indicates that the voltage is Rayleigh distributed and the slope indicates the value 
of the Rayleigh parameter, a 2 ~ 1 0.3 V 2 . 


CONCLUSIONS 

The system discussed in the previous sections for simulating the effects of space-to-space and 
space-to-ground optical communication channels is relatively easy to construct and enables 
controlled tests to be performed in the laboratory. The ultimate utility of, the simulator is 
dependent on the validity of the theoretical models chosen to represent each channel. When 
coupled with a laser source, the Channel Simulator is capable of providing l.og normal, Ray- 
leigh, Ricean, chi-squared or beta irradiance modulation and voltages as well as Rayleigh and 
Ricean modulated electric fields. It is intended that the simulator discussed in this paper will 
be incorporated into a high data rate optical communication system for laboratory testing of 
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Figure 8. Statistical test of beta channel. 
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the composite system performance. Comparisons with computer generated analyses will then 
be made to evaluate the accuracy of current analytical models for laser communication links 
undergoing fading. 
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